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We calculate the classical gluon field created at early times in collisions of large nuclei at high 
energies. We find that the field is dominated by the longitudinal chromoelectric and chromomagnetic 
components. We estimate the initial energy density of this gluon field to be approximately 260 
GeV/fm^ at RHIC. 
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Experiments are being carried out at the Relativistic 
Heavy Ion Collider (RHIC) and soon will be at the CERN 
Large Hadron Collider (LHC) to create and study quark 
gluon plasma. Data from RHIC indicate that in colli- 
sions of gold nuclei at ^/SNN — 200 GeV energy densities 
far in excess of the critical value required for deconfine- 
ment (cc « 2 GeV/fm'^) are reached [l|. Furthermore, 
the partonic phase seems to be thermalized after a very 
short time tq < 1 fm/c. While the evolution of the quark 
gluon plasnra in equilibrium can be described by rela- 
tivistic hydrodynamics [2]; the initial soft interactions of 
the nuclei and the thermalization process before the time 
To are still not completely understood. 

It has beeir argued that the initial dynamics for the 
collision of two very high energy nuclei is determined by a 
universal phase called the color glass condensate (CGC). 
This idea is based on gluon saturation at a scale Qs ^\d, 
la, ly, |7| . Slowly evolving and randomly distributed color 
charges in the nuclei are the sources of this gluon field. 
A simple implementation is the McLerran-Venugopalan 
(MV) model 3, 4] in which the gluon field is given by the 
solution of the classical Yang- Mills equations. 

In this Letter we calculate the gluon field at early times 
after the collision in the framework of the McLerran- 
Venugopalan model. We use an expansion of the Yang- 
Mills equations in powers of the proper time r. This is a 
near-field approximation which may be the most appro- 
priate use of the color glass coirdensate picture. We also 
estimate the initial energy density at the time of overlap 
of the nuclei using a simple model for the nuclear gluon 
distribution and coarse-graining methods to avoid ultra- 
violet (UV) singularities. More details and a discussion 
of applications will be provided elsewhere |8t] . 

In high energy collisions the two colliding nuclei are 
highly Lorentz contracted; therefore, the valence and 
large- a; partons are described by infinitcsimally thin 
sheets propagating on the light cone. Although each nu- 
cleus is color neutral as a whole, local color fluctuations 
do occur. At the moment of overlap, the color distribu- 
tions in nucleus 1 (-1- light cone) and 2 (— light cone) 
are pi(xj^) and /02(xj_), respectively. We use light cone 
coordinates x^ = (x" ± x^) /\/2 and x^ ::= {x^,x^). The 



distributions pk = p%t°' (fc = 1, 2) are functions with val- 
ues in 5[/(3). Since they resemble fluctuations of color 
we have to take the ensemble average of all allowed func- 
tions Pk at the end. It is convenient to choose an ax- 
ial gauge defined by A'^x^ + A~x'^ — 0. In this gauge 
the current generated by the charges pk takes the form 
J* = 5{x^)gpi^2{^±) and J' = 0, satisfying the equation 
of continuity [Z?^, J^] = 0. The gluon field generated by 
this current can be obtained by solving the Yang-Mills 
equations [D^,F^'•'] = J". 

The gauge potential A^ is a smooth function of x^ 
except for lines with propagating charge. We follow the 
authors of ref. Ijl who showed that the ansatz 
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{x-)e{-x+)a\ix^) + 9{x+)e{-x-)ati{xi_) 



+ e{x+)e{x-)ai{T,xi_). 
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satisfies the Yang-Mills equations in the different sectors 
of Minkowski space. Here upper Latin indices i, j, ... 
always refer to transverse components. The a\ and a^ 
are the purely transverse gauge potentials of nucleus 1 
and 2, respectively. They can be written with the help 
of transformation matrices Uk — e''^'=, k — 1, 2, such that 
they are gauge transformations of the vacuum: 



aUx±) = -U^^d'Uk,with 
Vi0fc(x_L) = 5^pfc(x^). 



(3) 
(4) 



a and a\ describe the field in the forward light cone 
(x+ > 0, x~ > 0) which is generated in the collision. 
They are smooth functions of xj^ and the proper time 
T = V2a;+a;^. They are independent of the space-time 
rapidity r/ = l/21n(a;+/x~) because the current J'^ is 
boost-invariant. In the forward light cone the Yang-Mills 
equations can be rephrased as [J]: 



drT^dra-[D\[D\a]] =0, 
D\dral]- igT[a,dra] =0, 



1 



drTdral - igT^[a, [D\a]] ~ [D\F^'' 



(5) 
(6) 
(7) 



a and a\ are connected with the single nucleus fields a\ ^ 
via boundary conditions at r = 0: 






(8) 
(9) 



A summation of the recursive solution in closed form does 
not seem feasible. However, we assert that an analysis 
using just the first few orders in r is extremely useful. 

The non- vanishing components of the field strength for 
the lowest three orders in r are 



An explicit analytic solution of ©-([T]) is not known. 
However, lowest order perturbative solutions [J| as well 
as numerical solutions [OJ, llOi are available. 

12| will even- 



F, 



(0) 



«5 ai,"2 



11 









2^/2 



([^(0)'-^m]=^ [-0(0) :^( 
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Decoherence and pair production 
tually destroy the classical field and lead to thermaliza- 
tion. Typical values for the thermalization time tq used 
in hydrodynamic calculations range from 0.15 fm/c to 1.0 
fm/c [lj,|j|. It is clear that the classical description breaks 
down before Tq. Hence, what we can hope to calculate in 
this particular framework is the short-term behavior of 
the gluon field, i.e. the near-field close to the light cone. 
The functions a and a\ are regular at r = 0. Therefore Here e*-' is the antisymmetric tensor. Thus the longitu- 
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Fi^^ - -[-Dfo)' [-0(01,-^(0)] 
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(15) 
(16) 

(17) 
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it is legitimate to solve the Yang-Mills equations using a 
power series in r. We write 



oc 

a(r, x_l) = ^ T"a(„) (x_l) 



(10) 



and similarly for a\. We also use expansions for the field 
strength tensor and covariant derivative in the forward 
light cone with coefficients F*f'\ and F)^^., respectively. 

Using these expansions in Eqs. ([5]) through (O yields 
an infinite set of equations for the coefficients «(„) and 
'-^\{nY ^^ lowest order in t, the fields are just given by 
the boundary conditions ([5]) and ([5]), 



^3(0) 



J=2 ; 



>^(o) = -»5K,a^]/2. 
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It is now possible to give a solution for any order in r 
recursively. It is straightforward to prove that for n > 1 



1 



*(n) 



n(n + 2) 



E 



D 



(k) 



[^(i):"( 



m)J 



fc+Z+m— n— 2 



\kM=n-2 



ni pi' 



(13) 



'W E 

k-\-l-\-rn—n—4: 



(fc), [-D(;),a(m 



One can immediately conclude that these fields vanish 
for all odd powers of r: a(2„+i) = 0, aLg^,]^-, = 0. Sim- 
ilar recursion relations hold for the field strength tensor. 
For brevity we only cite the relation for the longitudinal 
chromoelectric field E, = _F"' which is 






^—' n(m + 2) 

k+l+m=n-2 ^ ' 
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dinal chromoelectric field Ez and the longitudinal chro- 
momagnetic field B^ = F"^^ start with finite values at 
r = 0. The transverse electric and magnetic fields, which 
are linear combinations of the components F , are zero 
at T = and start at order r^. Generally, longitudi- 
nal fields have only contributions from even powers in r, 
transverse fields only consist of odd powers in r. 

This observation leads to the following space-time pic- 
ture. Inside the nuclei the color sources create purely 
transverse fields F'* = 5{x^)a\ 2 on the light cone. This 
is completely analogous to the abelian case. After nuclear 
overlap, non-abelian interactions between these fields cre- 
ate strong longitudinal chromoelectric and chromomag- 
netic fields, while the onset of transverse fields in the 
forward light cone is delayed. The situation resembles 
a capacitor with a longitudinal field, but it is important 
to realize that only the non-abelian nature of the gluon 
field can generate such a field for recoilless charges reced- 
ing from each other with the speed of light. The strong 
longitudinal fields at early times are an immediate conse- 
quence of the equations of motion; however, this fact has 
not received much attention before. Recently the strong 
pulse of longitudinal fields and its possible consequences 
have been discussed l^, [ij, llJ, lla] ■ 

In the second part of this Letter we would like to use 
our results to discuss the initial energy density eq = {T?^-^ ) 
for r ^ 0. Here T^'^ is the energy momentum tensor of 
the classical field which we expand in powers of r as well. 
We postpone all further discussion to a later publication 
[sl . Note that the recursion formulas use the gluon fields 
a\ as the starting point. However those have to be de- 
termined by solving the Yang-Mills equations dS]) , (|4]) for 
a single nucleus which is a difficult task. For our discus- 
sion here, difficulties associated with non-linearities in 
the boundary conditions are simplified by a mean-field 
approximation. As we will argue below it still represents 
the essential physics of the full solution. We achieve this 



by replacing Eqs. ([3]) and (|4]) with 



V^ - — 



= g^Pk. a'fc-— a>fc. (20) 
9 



The solution is formally linear in p^. However, we intro- 
duced a screening length R^ as a parameter which will 
depend on the charge distribution pk- 

The idea behind this approximation to Eqs. ([3]) and 
(HI is as follows. To lowest order in the charge density 
Pk we have a\. — — (l/.g)9Vfe + ^(Pk)- ^^'^ primary 
effect of the non-linearities is a partial screening of the 
field on length scales ~ l/Qs [5]. However, the screening 
is incomplete unless confinement is enforced in addition 
[l6|. In our approximation screening is provided by the 
scale Re ~ l/Qs and the result is perfectly infrared safe. 
For a reasonable estimate of the screening effect one can 
invoke the analogy to electric screening in a QCD plasma 
which implies R^^ — AusG /i where a is the area density 
of the number of color charges. 

Let us now consider the two nuclei as being made up of 
two ensembles of discrete charges, given by SU{i) matri- 
ces Ti_„ and T2^u at transverse positions b„. The charge 
densities can be written as 



Pfe(x±) 






i?(|xj 



bJ)T, 



k.u • 



(21) 



The index u represents coarse-grained cells in both nu- 
clei with a certain number of color sources, iVi.u and 
N2.U, respectively, in each of them, and R{r) is the spa- 
tial distribution of the effective color charge in each cell 
normalized to one. Note that the factorization of spatial 
and color degrees of freedom makes sense if the coarse- 
grained cells have sizes much smaller than 1/Qs- 

Then from Eqs. ([3]), (|4]) it follows that the gauge po- 
tentials can be written as 



aU^±)=9Yl 



X_L 



b„ 



■ffe,„G(|xi-bJ) (22) 



where G is the linear field of a single charge, i.e. G{r) — 
dip/dr with V5_(^(r) = R{r). The fk,u are SU{3) valued 
functions and can be interpreted as the modifications of 
the charges T^ ^ through non-linear interactions. One 
can expand Tk,u{r) — Tfe.u + 0{T^), where the first term 
is the linear contribution, and higher order terms reflect 
the screening from interactions with charges in neighbor- 
ing cells. We can now simplify the situation by applying 
the approximation ([20| . It amounts to the replacement 
Tk,uG -^ Tk^uG. The screening present in T is now fac- 
torized into a modified field profile G with 



vi 



i?2 



(p{r) = R{r) , G{r) = dip/dr . (23) 



In addition, we have to impose an ultraviolet (UV) cut- 
off Qq. The McLerran-Venugopalan model does not pro- 
vide a UV finite answer for the energy density at r = 0, 



as also noticed in 17|. This comes from the fact that 
hard modes with momentum much larger than Qs are 
treated correctly. They are better described by hard per- 
turbative processes. Therefore Qq should be chosen to be 
the cutoff between hard processes involving modes with 
transverse momentum pt > Qo, and the bulk modes with 
Pt < Qo- We realize that our coarse-graining provides 
this cutoff. To be more precise we choose Gaussian pro- 
files R for each charge with a width A = 1/Qo- Putting 
everything together we find the modified field profile is 
approximately 



Gix^) 



1 — exp(— x^/A^) 

2TTRr 



^^'1^ 



(24) 



NL 



NL 



and N^ ^ are the num- 



where Ki is a modified Bessel function. 

To calculate the expectation values of observables we 
discretize the functional integrals over the charge dis- 
tributions Pk and replace them with integrals over the 
group SU{3) at each point b„. The correct weight func- 
tion to be used for the integral for cell u in nucleus 
k is WN,jTk,u) = {Nj{TTNk,u)y expi-N,TlJNk,u) 
where iVc = 3 [l8|. As a straightforward generaliza- 
tion of the results in [l8| we use Nk^-, 

Ca/CfNI^ where iV,^^„, N^ 
ber of quarks, antiquarks and gluons in each cell. We 
define the area charge density in nucleus k as ak,u — 
^fc,u/(area of the cell u) for each cell. It is then straight- 
forward to define a continuous charge density (Tfe(xj^). 
The saturation scale usually contains the strong coupling 
and we set Q^ — as<J. 

To summarize, our model for the gluon field of a sin- 
gle nucleus deviates in two ways from the McLerran- 
Venugopalan model. First, for simplicity we use a 
mean-field approximation which reproduces the essential 
physics. Second, we implement a UV regulator Qo- To 
compare with the existing literature one can compute 
some quantities in the limit Qo ~^ 00 • For the correla- 
tion function of charges one obtains (p^(x^)p^(yj^)) — 
'^"^(o'fc/S) (5-^(xj^ — y±) for constant densities (Tk (here 
a, 6 = 1, . . . , iV^ - 1) 0, i, i, [l3|. In the same limit we 
find that the field correlator {A'^{x±)A'^{y±)) is a good 
approximation of the analytic form [5| . 

We are now ready to use our model of the single nu- 
cleus gluon field to obtain the initial electric and mag- 
netic energy densities e^ = (Tr(Fir)2) and em = 
{Ti{F?^-.)^). After evaluating the expectation values 

(z2Tr([ri,„,r2,„][Ti,„,,T2,,-])) = Suu'Svv'Ni,uN2.,/Ne 
(the trace refers to color) we find 



eij(x±) = |r XI ^i,«^2,^ ((x± - b„) • (x^ - b„)) 
G(|xi-b„|)2G(|x^-b„|)2 



-bJ2 



|x± 



bJ2 



(25) 



For tM the square on the first Hne has to be replaced by 
(e*^ {KK - x±{K - bi))Y. We evaluate this result for 
the center (xj^ = 0) of two large nuclei with radius Ra 
colliding head-on, so that Ra S> i?c S> A. The contri- 
butions of the two nuclei to Eq. ([25|) factorize if the size 
of each cell is small, as was also noticed in [l7| . Assum- 
ing that the charge densities ak are roughly constant in 
the center of each nucleus we find to good approximation 
that 



££ 



en - 2^0 



^a,a,\n'{l + ce). (26) 



This result only depends on the ratio of scales C = Rc/^ 
and c « 0.42 is a numerical constant. 

The charge densities <Tk , whose fluctuations create the 
color distributions pk, are given by the large- a; partons 
in the nuclei. To give a numerical estimate for two Au 
nuclei colliding at RHIC energy we count all partons in 
the nuclei above the cutoff scale Qo, similar to the pro- 
cedure in [19|. In practice we determine a = ai — a2 a,s 
a function of Qo using CTEQ parton distributions. Note 
that while the screening length Re in a nucleus is a phys- 
ical quantity, the cutoff Qo is unphysical. We observe 
that Qs = y/asCT ~ 1/Rc is indeed almost independent 
of Qo; however, the energy density eq is not. The resid- 
ual logarithmic dependence on Qo should vanish if we 
match classical and hard perturbative results in the re- 
gion where they are comparable 19]. 

Fig. [1] shows our estimate for the inital energy den- 
sity eo in the center of the collision as a function of 
the UV cutoff Qo for central collisions at RHIC us- 
ing ("^ = Ql/{as(j). We find that Q^ only varies be- 
tween 1.4 and 1.7 GeV if Qo is varied between 1 and 
10 GeV. For a reasonable cutoff Qo = 2.5 GeV we have 
eo ~ 260 GeV/fm'^. This is compatible with a value of 
130 GeV/fm3 at r = 0.1 fm/c found by T. Lappi in [Hi- 
Note that eo only takes into account gluon modes with 
transverse momentum less than Qo. Results for finite r 
and the matching with hard processes to compute the 
total energy density and to eliminate the sensitivity to 
Qo will be discussed in a forthcoming publication [8J . 

To conclude, we introduced a near-field expansion to 
solve the classical Yang-Mills equations for the collision of 
two large nuclei in the color glass picture. We found that 
strong longitudinal chromoelectric and magnetic fields 
dominate at early times. Using a coarse-graining of color 
charges we derived a simple expression for the initial en- 
ergy density of the soft gluon field. A rough estimate 
implies values of about 260 GeV/fm'^ at r = for the 
center of two colliding nuclei at RHIC. 
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FIG. 1: Initial energy density eo for t — > and saturation 
scale Qs at RHIC as a function of the UV cutoff Qo. 
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